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ABSTRACT 

We present the results of a joint observational campaign between the Green Bank radio telescope 
and the VERITAS gamma-ray telescope, which searched for a correlation between the emission of 
very-high-cncrgy (VHE) gamma rays {E^ > 150 GeV) and Giant Radio Pulses (GRPs) from the 
Crab pulsar at 8.9 GHz. A total of 15366 GRPs were recorded during 11.6 hours of simultaneous 
observations, which were made across four nights in December 2008 and in November and December 
2009. We searched for an enhancement of the pulsed gamma-ray emission within time windows placed 
around the arrival time of the GRP events. In total, 8 different time windows with durations ranging 
from 0.033 ms to 72 s were positioned at three different locations relative to the GRP to search for 
enhanced gamma-ray emission which lagged, led, or was concurrent with, the GRP event. Further, 
we performed separate searches on main pulse GRPs and interpulse GRPs and on the most energetic 
GRPs in our data sample. No significant enhancement of pulsed VHE emission was found in any of 
the preformed searches. We set upper limits of 5-10 times the average VHE flux of the Crab pulsar on 
the flux simultaneous with interpulse GRPs on singlc-rotation-period time scales. On ~8-second time 
scales around interpulse GRPs, we set an upper limit of 2-3 times the average VHE flux. Within the 
framework of recent models for pulsed VHE emission from the Crab pulsar, the expected VHE-GRP 
emission correlations are below the derived limits. 
Subject headings: gamma rays: stars, pulsars: individual: B0531+21 
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1. INTRODUCTION 

The Crab pulsar, PSR B0531+21, is a powerful young 
pulsar and one of the most studied objects in the sky. 
It is one of the b rightest pulsars in the high-energy 
gamma-ray regime (jFierro et al.|[l99l lAbdo et al.|[2010l) 
and the onl y puls ar so far to be detected above 100 GeV 
(jAliu et alj(2011l : lAleksic et al.ll2012[ ). The Crab pulsar 
is also one of only several known pulsars which exhibit 
the giant radio pulse phenomenon ({Knight 2 00^ : single 
radio pulses with flux densities that greatly exceed the 
average pulse flux density and can, at their maximurn , 
be as bright as a few million Janskys ()Soglasnovll2007[ ). 
The energies of giant pulses follow a power-law distribu- 
tion (|Cordes et al. 2004; Popov & Stappcrs 2007.) m con- 
trast to the regular pulses that obey a Gaussia n or log- 
normal distribution (jBurke-Spolaor et al.l 120121 ). High- 
resolution time measurements of individual GRPs reveal 
pulses whi ch can be a few micro seconds to a few nanosec- 
onds wide (jHankins et al.l2003D , with the narrower GRPs 
possessing the highest flux density. 

The pulse profile of the Crab pulsar is dominated 
across the electromagnetic spectrum, from radio waves 
to VHE gamma rays, by two emission peaks, referred 
to as the main pulse (MP) and interpulse (IP). These 
peaks occur at phases 0.0 and 0.4, respectively, with 
a "bridge" of enhanced emission appearing between 
these two peaks in the optical, X-ray and gamma- 
ray pulse profiles (lOosterbroek et al.ll2008t iMineo et al.l 
[20061: IKuiper et al.l 120011: lAbdo et al.l I201(E Above 
100 GeV, significant emi ssion is oiily observed during 
the m ain and inter-pulse ()Aliu et al.ll2011l : lAleksic et al.l 
l2012f l. In the radio-pulse profile several different compo- 
nents appear at different frequencies, including a precur- 
sor to the main radio pulse and two "high frequency com- 
ponents" (HFCs), which only appear above ~5 GHz and 
occur at phases '^0.77 and '^0.93 ([Moffett fc Hankini 
|1996[) . In the Crab pulsar, giant radio pulses have been 
observed to occur during both the main pulse and the 
interpulse. Indeed, it has been suggested that the Crab 
radio profile is composed of entirely GRPs, with the 
"regular" pulsar pulse corresponding to the main pulse 
precursor component ijPopov et al.ll200"6l ). The case for 
the regular pulse is difficult to disentangle, since sin- 
gle regular pulses cannot be observed over the nebular 
background. Stark differences in the characteristics of 
main pulse GRPs and interpulse GRPs have been re- 
vealed by observations of the Crab pulsar between 6 and 
10.5 GHz. In this frequency range, interpulse GRPs are 
typically several microseconds long and populat e a set 
of re gularly spaced frequency bands (jHankins fc Eilekl 
|2007| ). On the other hand, main pulse GRPs exhibit 
broadband spectra and appear as a succession of nar- 
row pulses ranging from unresolve d widths below 0.4 ns 
to widths of a few m icroseconds (jHankins et al.l 120031 : 
iHankins fc Eilekl 120071) . These striking differences sug- 

Institute of Technology, Bishopstown, Cork, Ireland 

30 ASTRON, The Netherlands Institute for Radio Astronomy, 

Postbus 2, 7990 AA, Dwingeloo, The Netherlands 

3^ Astro Space Center of the Lebedev Physical Institute, 

Profsoyuznaya str. 84/32, Moscow 117997, Russia 

3^ To whom correspondence should be addressed. E-mail: 

mccann@kicp.uchicago.edu (A.M.); nepomuk.otte@gmail.com 

(A.N.O.); schroedter@veritas.sao.arizona.edu (M.S.) 



gest that the emission mechanisms may differ between 
the main pulse GRPs and interpulse GRPs in the Crab 
pulsar above 5 GHz. 

The mechanisms responsible for the generation of 
GRPs are still unknown. Changes in the coherence of the 
plasma beam, which is believed to be responsible for the 
normal pulsed radio emission, can in principle explain 
the generation of GRPs. However, mechanisms which 
primarily affect the coherence of the plasma have no ef- 
fect on the incoherent emission from pulsars and, thus, 
no enhancement is expected in the gamma-ray emission 
in connection with GRP events. Mechanisms which in- 
crease the rate of particle production within the magnc- 
tospheric emission region, or which change the direction 
of the emission beam, should, however, affect the higher 
energy incoherent emission. Such mechanisms may cre- 
ate an enhancement in the gamma-ray emission from the 
pulsar. 

The exceptional emission band structure seen in inter- 
pulse GR Ps from the Crab pulsar at high frequencies re- 
ported bv IHankins &: EileH (j2007[ ) prompted a qu antita- 
tive model for their generation bv iLvutikovl (j2007[ ). This 
model argues that GRPs are generated in a dense plasma 
region close to the last closed magnetic field line. Oc- 
casional magnetic field reconnections excite the plasma 
producing a particle beam with a high Lorentz factor 
which emits a GRP via anomalous cyclotron resonance. 
This model, when fit t o the interpulse GRP data from 
IHankins fc Eiiei3 (j2007l ). returns a predicted Lorentz fac- 
tor for the beam, 7 = 0(10^), which is large enough to 
cause the particles within the beam to generate curvature 
photons with energies as high as tens of GeV. Thus, a 
feature of this model is a prediction of an enhancement in 
the gamma-ray emission correlated with high frequency 
interpulse GRPs from the Crab pulsar. 

Several studies have been performed to examine a pos- 
sible connection between G RPs in the Crab and h igher- 
energy incoherent emission. iLundgren et al.l (jl995[ ) stud- 
ied 50-220 keV gamma rays recorded by the OSSE instru- 
ment on board NASA's Compton Gamma-ray Observa- 
tory along side Green Bank radio data recorded at 800, 
812.5 and 1330 MHz. No enhancement was seen in the 
gamma-ray emission on 2-minute time scales concurrent 
with GRP events. This study yielded an upper limit 
on the average gamma-ray flux from the pulsar over 2 
minutes concurrent with GRPs of 2.5 times the average 
pulsed gamma-ray flux. Based on the non-detection of 
a gamma-ray enhancement, the authors argue in favour 
of coherence changes within the emission plasma as the 
source of GRP generation. S hearer et al.l (j2003[ ). how- 
ever, observed a signiflcant 3% enhancement in the op- 
tical main pulse concurrent with the period of emission 
of GRPs measured at 1380 MHz. No enhancement was 
seen in the interpulse. This observation suggests a link 
between coherent radio emission and incoherent optical 
emission in the Crab pulsar. Small changes in the pair- 
creation rate leading to localised density increases within 
the emission plasma could create the GRP event and 
provide a small enhancement in the optical incoherent 
emission. 

Recentlv lBiTous et al.l (j20Tl studied 0.1-5 GeV gamma 
rays recorded by the Fermi-hAT in conjunction with 8.9 
GHz observations of the Crab performed by the Green 
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Bank radio telescope. No enhancement was seen in the 
gamma-ray flux within the main pulse, interpulse and 
bridge, during single rotation periods which contained 
a GRP. This yielded a limit on the flux from the pul- 
sar during rotation periods concurrent with GRPs of less 
than 4 times the average gamma-ray flux from the pulsar. 
This result suggests that enhanced pair-creation is not 
the dominant factor in the emission of GRPs in the Crab 
and supports the idea that local coherence changes in the 
magnetospheric plasma cause GRPs, or, that the puta- 
tive enhanced gamma-ray emission is beamed in a differ- 
ent direction than the radio emission and th us unobserv- 
able f rom Earth. Another recent study by iBilous et al.l 
()2012D . which examined 1.4-4.5 keV X-rays from Chan- 
dra observations of the Crab pulsar in conjunction with 
radio observations by the Green Bank telescope at 1.1-1.9 
GHz, found no enhancement in the X-ray flux in connec- 
tion with GRPs. This study yielded limits on the flux 
enhancement in the main pulse and interpulse concurrent 
with GRP events of 10% and 30%, respectively. 

Along with the aforementioned studies which searched 
for a correlation between GRPs and incoherent emission 
below 10 GeV, a GRP correlation search was performed 
on the Crab pulsar at TeV energies in the early 1970s 
(jArgyle et al.l 1197^ . This pioneering attempt was per- 
formed with the 25-m DRAG radio telescope and the 
Whipple gamma-ray telescope almost two decades be- 
fore the first conclusive det ection of a TeV gamma-ray 
source (jWeekes et al.lll989( ). No TcV-GRP correlation 
was observed and we are aware of no other study on this 
subject with subsequent TeV instruments. 

VERITAS is an array of four imaging atmospheric 
Cherenkov gamma-ray telescopes located at the Fred 
Lawrence Whipple Observatory a t the base of Moun t 
Hopkins in southern Arizona (i Holder et all [20061 ). 
Each telescope u s es a 12 m Davics-Cotton reflector 
(jPavies fc Cotton! Il957t) instrumented with a photo- 
multiplier tube (PMT) camera with 499 pixels. The 
array is sensitive to gamma rays with energies between 
100 GeV and 30 TeV. The recent detection by VERI- 
TAS of pulsed em ission above 100 GeV from the Crab 
(jAliu et al.ll2011| ) is not expected within the context of 
contemporary pulsar models in which curvature radia- 
tion is the dominant gamma-ray emission mechanism. 
The origin of the power-law extension of the spectral en- 
ergy distribution above 100 GeV seen now by both VER- 
ITAS and MAGIC (Alcksic et al. 20l2) is poorly under- 
stood. A possible explanation is that the VHE emis- 
sion is caused by inverse- Compton (IC) up-scattering of 
soft photo ns by pa.rticles accelerated in the outer magne- 
tosph ere (jRomanil 119961: iLvutikov et~all 120121: iDu et al.l 
l2012r ). Another possible explanation is that the pulsed 
VHE emission is produced when pulsed magnetospheric 
X-ray photons are up-scattered by re lativistic plasma 
in th e wind outside the light cylinder (jAharonian et al.l 
120121) . Given that the origin of the VHE emission and the 
generation of GRPs in the Crab pulsar are both poorly 
understood energetic phenomena, we are motivated to 
probe whether they are connected. 

The remainder of this paper is structured in the follow- 
ing way. In §2 we describe the simultaneous observations 
of the Crab pulsar made with the GBT and VERITAS 
gamma-ray telescopes. We also describe the extraction 
of GRP events from the radio dataset. In §3 we dis- 
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Pulsar Phase 

Figure 1. Summary of the gamma-ray and radio observations. 
Panel (a) shows the gamma-ray pulse profile for the pu blished 107 
hour VERITAS Crab Pulsar dataset HAliu et al.|[20nD . The over- 
laid solid line is the best fit function determined from a maximum 
likelihood fit to the unbinned VERITAS data. Panel (b) shows the 
gamma-ray (grey) and giant radio pulse (hatching) profiles for the 
four nights of simultaneous observations. Panels (c),(d),(o) and 
(f) show the arrival time at the barycenter versus phase for the 
observation dates 54829, 55153, 55158 and 55180 in MJD respec- 
tively. The grey-scale histograms show the VERITAS data while 
the blue, green and red points are the radio data with energies less 
than 100 Jy/^s, greater than 100 Jyfis and greater than 150 Jy/is, 
respectively (see colour figure in on-line version). 

CUSS the gamma-ray data processing and our strategy for 
searching for correlated emission within the radio and 
gamma-ray datasets. In §4 we describe a Monte-Carlo 
simulation of the Crab pulsar signal within the VERI- 
TAS dataset, which is used to calculate the strength of 
any correlation found in the datasets. In §5 we detail the 
results of the enhancement search and in §6 we provide 
some concluding remarks. 

2. OBSERVATIONS 

Simultaneous observations of the Crab pulsar were 
made by VERITAS and the 100-m Robert C. Byrd Green 
Bank Telescope (GBT) on the 29th of December 2008, 
the 18th and 23rd of November 2009 and the 15th of De- 
cember 2009. A total of 11.6 hours of simultaneous data 
was recorded across these four nights with 2.3, 3.33, 3.04 
and 3.00 hours acquired, respectively. Table [1] gives a 
summary of the data used in this study. 
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Date # GRPs # GRPs # 7-ray T Overlap Center Freq. A samp. 



[MJD] 




Overlap 


Candidates 


[min] 


[MHz] 




54829 


280 


150 


5566 


138.365 


8832 


40.96 


55153 


7813 


5937 


8738 


200.232 


8900 


3.2 


55158 


3771 


3243 


7928 


182.881 


8900 


3.2 


55180 


6916 


6036 


7861 


180.107 


8900 


3.2 


Total 


18780 


15366 


30093 


701.585 







Table 1 

A summary of the radio and gamma-ray data used in this study. Measurements were made on four separate nights with a total 11.6 hours 
of simultaneous data accumulated. The significantly smaller number of GRP events found on the first night of observations is likely due 

to higher attenuation of the radio input signal. 



2.1. GBT Observations and GRP selection 

The radio data presented here were acquired using 
the Green Bank Ultimate Pulsar Processor Instrument 
(GUPPI) in search mode. The total bandwidth of 
800 MHz was centered at 8832 or 8900 MHz and split 
into 128 or 256 frequency channels for our observations 
in 2008 and 2009, respectively. Full Stokes parameters 
were recorded at a sampling interval of 40.96 /is in 2008, 
whereas total intensity was recorded with a sampling in- 
terval of 3.2 us in our 2009 observations. 

Recorded data from every session were dedispersed 
using the PRESTO pulsar software packageo, and 
searched for all single-pulse events with a peak signal- 
to-noise ratio seven times greater than the average radio 
signal. The dedispersion was done using the contempo- 
raneous dispersion measure (DM) value stated in the Jo- 
drell Bank Crab pulsar monthly ephemeri^3 (|Lvne et al.l 
119931. The GRP selection was performed with the 
PRESTO tool singlepulse_search.py, which convolves 
the dedispersed time series with a series of boxcar func- 
tions of different widths. The times of arrival of the se- 
lected GRP events were converted into TempcP^ format 
and transformed to barycentric dynamical time (TDB) 
for the correlation analysis with the VERITAS events. 

Between the observation sessions in 2008 and 2009, the 
DM has increased from 56.7883 pc cm^'^ in December 15, 
2008 to 56.8279 pccm-^ in December 2009. We estimate 
timing errors at our high observing frequency of 8.9 GHz 
to be ~ 0.4 /xs, which is less than our sampling time. 
Here we assume a smooth variation in the DM and that 
the changes in the DM between our observing sessions 
are certainly less than the measured change of the DM 
of « 0.04 pccm~^ over the course of one year. 

The system equivalent flux density is mostly deter- 
mined by the Crab Nebula. Flux densities of the Crab 
Nebula were calculated with the relation S{f) = 955 x 
(//GHz)-°-27 Jy (ICordes et al.l[200l . accounting for the 
fact that at 8.9 GHz the solid angle of the GBT beam 
covers only 6.25% of the area occupied by the nebula. We 
estimate the system equivalent flux density to be about 
0.2 Jy for our 2008 observations, and about 0.7 Jy for the 
observing session in 2009. The smaller number of GRPs 
found in the 2008 dataset is likely due to a higher atten- 
uation in the receiver system in 2008 where two back- 
end readouts were used together; the GBT Pulsar Spigot 
Card and the GBT Mk5 disk readout. The clumping of 
the barycentric arrival times of the GRPs clearly seen in 

http: //www. cv.nrao ■ edu/$\sim$sransom/presto/| 
^* http: //www. jb .man. ac .uk/p ulsar/crab .html| 
http: //tempo . sourcef orge .net | 



the 2009 observations (see Figures 1 and 4) agrees with 
what one would expect from refractive interstellar scintil- 
lations (RISS). At 8.9 GHz the characteristic time scale 
of RISS is about 80 minutes (Bilous et al. 201^. 

2.2. VERITAS Observations 

The VERITAS observations were made under the best 
possible sky conditions with each of the four telescopes 
fully operational. These observations were made with 
zenith angles smaller than 30 and an average zenith 
angle of 16 , yielding the lowest possible e nergy thresh- 
old. A ll data were acquired in wobble mode (jFomin et al.l 
I1994D . This mode of observation places the object under 

study at an offset of 0.5 from the centre of the field of 
view of the telescope, allowing for a simultaneous mea- 
surement of the background in other regions of the field 
of view with the same acceptance as the region contain- 
ing the source. 

Before an event is recorded with VERITAS, at least 
two of the four telescopes must trigger on a Cherenkov 
flash within 50 ns. A single telescope trigger is formed 
when three or more adjacent PMTs in the camera reg- 
ister at least six photoelectrons within 9 ns. When the 
trigger condition is satisfied, the PMT traces in each tele- 
scope are read out by 500 mega-sample per second flash 
analogue-to-digital converters which are located at each 
telescope. For each event, the VERITAS data stream 
contains the GPS time stamp from each of the four tele- 
scopes along with the digitised PMT-traccs. In the entire 
VERITAS dataset presented in this work, the four GPS 
timestamps never diverged by more than 10 /iS, providing 
a sound basis for correlation and timing analyses. 

3. DATA ANALYSIS 

3.1. Gamma-ray Event Selection 

Prior to performing a correlation search between the 
gamma-ray and radio datasets, the VERITAS data were 
passed through an analysis pipeline, which reconstructs 
the arrival direction and the energy of the gamma-ray 
candidate from the Cherenkov images recorded by the 
telescopes. For each event, the signal in each PMT is 
corrected for gain differences between the PMTs and 
signals which only contain noise are removed. Follow- 
ing this image cleaning, images which have a combined 
PMT signal (size) above 20 photoelectrons have their 
second moments calculated about their major and mi- 
nor axes ([Sillas 1985|). The arrival direction and impact 
parameter of the gamma-ray candidates are then calcu- 
lated from the intersection points of the major axes of 
the images from multiple telescopes when project ed on 
the sky and ground planes respectively ([Hofmann et alj 
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Figure 2. Positioning of the gamma-ray search windows with re- 
speet to the time of arrival of the GRP event at the solar system 
barycentre (SSB). For each search window duration, a window was 
positioned centered on the GRP, with leading search and lagging 
search windows placed directly before and directly after the cen- 
tered search window, respectively. 

11999( 1. Background suppression, including cosmic-ray re- 
jection, is performed by comparing measured event pa- 
rameters with Monte-Carlo gamma-ray simulations. Se- 
lection parameters include image brightness and shape, 
arrival direction, height of shower maximum, and impact 
distance, and are combine d in multi-dimensional en ergy- 
dependent look-up tables (jKrawczvnski et aT1l2006l ) . The 
optimal cut values were chosen a priori by modelling the 
gamma-ray signal from the Crab pulsar as a power-law 
with an integral flux of 1% of the Crab Nebula flux above 
100 GcV and a spectral index of —4. A grid search for the 
optimal gamma-ray cut values was then performed using 
this simulated signal with real Crab Nebula data as back- 
ground. The gamma-ray selection parameters and the re- 
sulting optimal cut values are: the angular separation be- 
tween the source location and the shower direction, theta 
(< 0.27 ), mean scaled width (<1.17), mean scaled length 
(<1.35), and height of shower maximum (>6.6 km). This 
cut optimisation yielded an energy threshold of 120 GeV 
for sources with a spectral index of -4. This gamma-ray 
analysis is identical to the analysis which was used in 
our paper on the first dete ction of the Crab pulsar above 
100 GeV (|Aliu et al.ll20Tll) . The GPS time of the candi- 
date gamma-ray events which passed cuts was converted 
to barycentric dynamical time and phase-folded using the 
Crab pulsar monthly timing ephemeris with an in-housc 
software package. The accuracy of these calculations was 
cross-checked using the Tempo program. 

3.2. Correlation Search Strategy 

The physical mechanisms which arc responsible for the 
emission of GRPs are not well understood. The nature of 
a possible connection between GRPs and the formation 
of high-energy incoherent emission from pulsars is also 
unknown. Given these unknowns, we prepared a corre- 
lation search strategy which probed for a connection be- 
tween GRPs and very-high-energy gamma-ray emission 
on different time scales, allowing for lagging or leading 
between gamma-ray and GRP emission. The searches 
were done by selecting only those gamma-ray events 
which arrived at the solar system barycenter within a 
certain time window around a GRP event. The dura- 
tions of the time windows were specified in units of pulsar 
rotations. We employed eight different time window du- 
rations lasting; 1, 3, 9, 27, 81, 243, 729 and 2187 pulsar 
rotations. This covers time intervals from 0.033 to 72.17 
seconds in time spacings which scale as log to the base 
3. Each time window was used to search for enhanced 
gamma-ray emission which lagged, led or was contem- 
poraneous with, the GRP. Thus, for each search window 
duration, a window was positioned centered on the GRP, 
with the leading and lagging windows placed directly be- 
fore and directly after the centered window, respectively 
(see Figure [2]). Finally, searches were performed consid- 
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Figure 3. The number of excess gamma-ray events in the phase 
interval -0.013 to 0.009 (main pulse) and 0.375 to 0.421 (inter- 
pulse) plotted against the total number of events for the complete 
107 hour VERITAS dataset ( Aliu et al. 20lJ). The excess growth 
is clearly linear. From a total excess of 1256±130 in 267088 events, 
there are 4.7±0.48 excess events in the VHE emission phases per 
1000 events selected. 

ering only main pulse GRPs, only interpulse GRPs and 
both main and interpulse GRPs combined. Thus, a total 
of 72 searches were performed. 

In each of the 72 searches, the phase of the gamma-ray 
events which fall within the search window is calculated. 
In our earlier mea surement of pulse d VHE emission from 
the Crab pulsar (jAliu et al.ll2011| ). we determined the 
phases of emission to be between —0.013 and 0.009 for 
the main pulse and between 0.375 and 0.421 for the inter- 
pulse. In this study we only consider VERITAS events 
that fall within these phase regions. When a gamma- 
ray event falls within the search window defined by a 
main pulse GRP, it has to lie within the VHE main pulse 
emission phase range for it to be considered in the en- 
hancement search. Corresponding selection criteria are 
applied to those gamma-ray events which fall within a 
window defined by an interpulse GRP. In the searches 
which consider both main and interpulse GRPs, gamma- 
ray events which fall within either of the VHE emission 
phase ranges are selected. The prescription for these 72 
searches was defined before the datasets were analysed. 

4. MONTE-CARLO TIME SERIES 

In order to determine the presence of an enhancement 
in the VHE emission from the pulsar correlated with 
GRPs, Monte-Carlo time series datasets were generated 
to model the gamma-ray data. Using a bin width of one 
second, raw trigger rate distributions (number of triggers 
per second) were compiled for each VERITAS observa- 
tion run. These distributions were used as probability 
density functions from which random event times were 
drawn and sorted, from earliest to latest, producing ran- 
dom time series with the same temporal characteristics 
as the real VERITAS data. Further, to enable the mea- 
surement of the level of any VHE flux increase seen in 
the VERITAS data, or to facilitate the calculation of flux 
upper limits, these Monte-Carlo time series datasets were 
injected with a simulated signal from the Crab pulsar as 
explained in the following section. 

From our earlier analysis of the complete 107 hour 
VERITAS Crab pulsar dataset, we determined that the 
excess of events, which fall within the VHE main and 
interpulse emission regions, grows line arly with r espect 
to the total number of events selected (jAliu et al.|[2011i) . 
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(a) Event Arrival Times (b) Event Separation 



Figure 4. Panel (a) shows the distribution of arrival times at the solar system barycenter of the gamma-ray (grey), radio (hatching) 
and Monte-Carlo data (circles) for the observations made on MJD 55153. The sharp drops in the gamma-ray data rate are due to gaps 
in observation between successive 20 minute exposures. Panel (b) shows the arrival time event-separation distribution for each pair of 
consecutive events in the complete gamma-ray and Monte-Carlo datasets. In both panel (a) and (b), the Monte-Carlo data contain an 
injected signal from the Crab pulsar at the level of the measured VHE pulsar flux {xi = 1) and are plotted as a profile histogram, with 
the height and error bar of each bin representing the mean and standard deviation determined from 500 simulated time series datasets 
respectively. 



Date [MJD] 


Total Number 
Nt 


Number of Excess 


Number of Background 


54829 


5566 


26.2 


5539.8 


55153 


8738 


41.1 


8696.9 


55158 


7928 


37.3 


7890.7 


55180 


7861 


37.0 


7824.0 



Total 30093 141.6 29951.4 



Table 2 

The number of signal and background events which are present in the VERITAS data recorded on the four nights of observation as 
determined by the linear relationship discussed in the text. The number of background events is defined as the total number of events 
minus the estimated excess. When generating a Monte-Carlo time series which models the observations made on a given date, the number 

of background events to be generated is drawn from a Poisson distribution with a mean Ni,g . The number of pulsar events to be 
generated is drawn from a Poisson distribution with a mean Nex x Xi, where Xi is used to scale the number of excess events to the desired 

pulsar flux level. 




0.1 0.2 0.3 

Pulsar Phase 

Figure 5. The VERITAS Crab pulsar phasogram determined 
from the 11.6 hours of simultaneous exposure together with pha- 
sograms determined from Monte-Carlo time series datasets. The 
Monte-Carlo data contain a simulated signal from the Crab pul- 
sar at various flux levels and are plotted as a profile histogram, 
with the height and error bar of each bin representing the mean 
and standard deviation determined from 500 simulated time series 
datasets respectively. 



This linear growth is plotted in Figure [3] With a total 
excess of 1256±130 in 267088 events we determine that 
there are 4.7±0.48 excess events in the VHE emission 
phases per 1000 events selected. Given this linear rela- 
tionship, the number of excess events which lie within 
the VHE emission phases, for any subset of the VERI- 
TAS data, can be estimated (see Table Further, we 
know that the emission peaks in the VERITAS Crab pul- 
sar phasogram can be modelled by two Gaussians sitting 
on a uniform background. This was determined by an 
unbin ned maxirnum li kelihood fit of the VERITAS phase 
data ()Aliu et al.ll2011t ). Using these two observations we 
can generate a Monte-Carlo time series inhabited by a 
simulated Crab pulsar signal at any desired flux level. 

If a given VERITAS observation has a total number of 
events, Nt, using the linear relationship discussed above, 
the number of excess pulsar events expected within this 
sample, N^x, can be estimated. The number of back- 
ground events is then Ni,g ~ Nt — N^x ■ Now assume we 
want to model this VERITAS observation with a sim- 
ulated signal from the Crab pulsar at a flux level Xi, 
where Xi is in units of the average pulsar flux level mea- 
sured with VERITAS. We draw a number of background 
events. N^g, from a Poisson distribution with a mean 
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Combined before 
Combined centered 
Combined after 
Interpylse before 
Interpylse centered 
Interpylse after 
Main pulse before 
t^ain pulse centered 
t^ain pulse after 
■ Range Expected from MC 



Searcfi Window Duration [Pulsar Periods] 

Figure 6. The number of gamma-ray events selected (symbols) 
versus the duration of the search window. The grey regions de- 
note 68% containment intervals about the mean of the distribution 
(black line) determined from searches performed on the Monte- 
Carlo data. The Monte-Carlo dataset is composed of 500 simulated 
time series containing an injected signal from the Crab pulsar at 
the level of the measured VHE pulsar flux {xi = 1). The absence 
of a symbol in 4 of the searches indicates that the number of se- 
lected events in each of these cases was zero. All searches return 
values which are consistent with what is found from searches on 
the Monte-Carlo datasets. The slight negative and positive shift of 
the x-position of the before and after symbols is done as a visual 
aid to prevent clutter about the common x-coordinate value for the 
before, centered and after symbols. 

Nhg and we draw a number of pulsar events, N^x, from 
a Poisson distribution with a mean N^^ x Xi. Now, us- 
ing the VERITAS raw data rate distributions, discussed 
earher, we draw Nt = N^g + N^x random arrival times. 
A fraction of these events, N^x/Nt, is randomly selected 
to contain the injected pulsar signal. This is done by 
shifting the arrival time of these events to the nearest 
time which, when barycentered and phase-folded, would 
correspond to a random phase value drawn from the 
double- Gaussian function which parametrises the VER- 
ITAS Crab pulsar phasogram. This small shift in time 
(<17 ms) has no effect on the overall data rate charac- 
teristics of the Monte-Carlo datasets, and is only applied 
to a small fraction of randomly selected events within 
each dataset. This procedure was used to generate sets 
of Monte-Carlo data, containing a signal from the Crab 
pulsar at a chosen flux level, which model all the data 
rate characteristics of the real VERITAS data, account- 
ing for the Poisson fluctuations inherent in the VERITAS 
measurement of the Crab pulsar flux. Examples of the 
match between the Montc-Carlo time-scries datasets and 
the VERITAS gamma-ray data are shown in Figures |4] 
and [3 

5. RESULTS 

5.1. Search for enhanced gamma-ray emission during 

GRPs 

Each enhancement search yields a number, iV, which 
is the number of VERITAS events that meet the specific 
search criteria, i.e., they fall within a time interval deter- 
mined by a GRP with a phase value within the required 
VHE-emission phase range. The Monte-Carlo time-series 
datasets were subject to the same search which was per- 
formed on the real VERITAS dataset. Given that we 
generated a large number of Monte-Carlo time series, 
a search on a Monte-Carlo dataset will yield a distri- 
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Figure 7. The deviation of number of events selected in the en- 
hancement searches on the real data compared to the mean number 
expected from the Monte-Carlo distributions in units of their stan- 
dard deviation. The positive shift of the mean of the distribution 
determined from searches on Monte-Carlo datasets with no pulsar 
signal flux signal injected (grey) is due to presence of the pulsar 
signal in the VERITAS dataset and the absence of a pulsar signal 
in the Monte-Carlo. Including in the Monte Carlo a Crab pulsar 
signal at the level of the average measured VHE pulsar flux yields a 
distribution (solid black), which has a mean compatible with zero 
and a standard deviation close to one. This indicates that there is 
no significant enhancement in the gamma-ray flux correlated with 
GRPs. We note that the searches are correlated because the same 
events can be selected in different searches. 

bution of the number of selected events, which is ap- 
proximately Gaussian, and will have a mean, ^, and a 
variance, ct^ = ^. The number of events selected in a 
given search, N, can be compared to the mean num- 
ber expected from the Monte Carlo when no enhance- 
ment is present, fj.. Such a comparison is plotted in Fig- 
ure [6l showing the number of events selected in searches 
of the VERITAS data, compared to the mean number 
selected from a Monte-Carlo dataset. Here, the Monte- 
Carlo dataset is composed of 500 simulated time series, 
each containing an injected signal from the Crab pulsar 
at the level of the average measured VHE pulsar flux 
{xi = 1). Using the formula, S = {N ~ fL)/a, one can de- 
termine the statistical significance of any deviation of the 
measured number of events, from the number expected 
in the absence of any enhancemenlF^. The significance 
distribution derived from the 72 searches is plotted in 
Figure [T] This distribution has a mean compatible with 
zero and a standard deviation close to one, indicating 
the absence of any significant enhancement in the VHE 
gamma-ray flux within the specified search windows po- 
sitioned around GRPs observed at 8.9 GHz. 

5.2. Calculation of the upper limit on the flux increase 

We follow the same prescription as iBilous et al.l (|2011[ ) 
to compute the upper limit on the VHE flux during the 
enhancement searches. From Bayes' theorem, the poste- 
rior probability that the pulsar flux is F given that we 

Twelve of the searches on the Monte-Carlo datasets yield 
distributions which are more Poissonian in shape than Gaussian 
(m ^ to), meaning that the strict equivalence of the significance 
formula breaks down. The resulting assertions are, however, un- 
changed. 
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Figure 8. Panels (a) to (h) show typical examples of the distribution of the number of selected events found in searches performed on 
main pulse GRPs. These distributions (grey) were determined from Monte-Carlo datasets containing an injected signal from the Crab 
pulsar at the level of the measured VHE pulsar flux. The solid black curve shows the best fit Poisson curve for each distribution. From 
the fits to distributions such as these, the values P(N\F = Xi) can be determined. Panel (i) shows the distribution of reduccd-x^ values 
(x^/ndf) for all of the Poisson fits to the Monte-Carlo distributions performed in this work, indicating that the Poissonian form accurately 
describes the data. 



observed N events in an enhancement search is 



P{F\N) 



P{N\F)P{F) 

OO 

/ P{N\F')P{F')dF' 





(1) 



where P{N\F) is the hkehhood of selecting A'' events in a 
search when the pulsar flux is F. For simplicity, all flux 
values specified here are cast in units of the average flux 
of the Crab pulsar measured with VERITAS. P{F), the 
prior distribution of the flux, is an uninformative prior 



which we set to be 



P{F) 



C if 1<=F <= 50 
if F < I,i^ > 50 



(2) 



where C is a non-zero constant. This means that during 
the emission of a GRP we consider the pulsar flux to be, 
with uniform probability, between I and 50 times the 
average Crab pulsar flux, and to have zero probability 
otherwise. 

From the searches performed on the Monte-Carlo 
datasets we can explore the likelihood value, P{N\F), for 
a range of flux values, by generating a likelihood curve. 
Each point in the curve is computed by probing the dis- 
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Figure 9. Panel (a) and (b) show the results of correlation searches centered on GRPs with search windows lasting 27 pulsar rotations. 
Panel (a) plots the posterior probability P(F = Xi\N) for a range of flux levels, given the number of events selected in the enhancement 
search. These probability distributions peak close to a flux value of 1, meaning that the flux is most likely unchanged during the GRP 
events. The error bars represent the uncertainty due to flnite statistics in the Monte-Carlo distributions and the fitting procedure used 
to determine the probability value from the these distributions. Panel (b) shows P{F < Xi\N')^ the cumulative distributions of the 
corresponding probability distributions shown in panel (a). From such curves, an upper limit on the level of the enhanced fiux can be read 
at a desired confidence level. In the given example, the 95% confidence level upper limits on the enhanced flux are 17.2, 3.5 and 2.9 times 
the flux measured by VERITAS for the main pulse, interpulse and combined searches respectively. 



tribution of the number of selected events determined 
from a Monte-Carlo dataset generated with a given flux 
level, F = Xi (see Figure [Sj for some example distribu- 
tions). The y- value of each point in the curve is the 
fraction of the Monte-Carlo datasets which, when sub- 
jected to the specific enhancement search, yielded the 
same number of coincident events as were found in the 
real data. In practice, we fit the distribution yielded 
from searches on the Monte-Carlo dataset with a Pois- 
son function, and determine the likelihood from the fitted 
function. This is done to minimise fluctuations caused 
by the finite statistics used to compile the Monte-Carlo 
distributions. Figure [5| shows some examples of the fits 
to the Monte-Carlo distributions along with the distribu- 
tion of the reduced-x^ values (x^/ndf) for all of the fits 
performed. From the reduced-x^ distribution, which has 
a mean value close to 1, it is clear that the Poisson func- 
tional form accurately describes the distributions derived 
from searches performed on the Monte-Carlo datasets. 

The likelihood curve for each search was compiled us- 
ing Monte-Carlo datasets with injected flux levels rang- 
ing from, = 1 to 50, with step sizes of 0.1 between 
1 and 10, 0.25 between 10 and 20, and 1 between 20 
and 50. Given our choice of the prior P{F), and that 
we evaluated the likelihood at an array of discrete flux 
levels, Equation 1 for the posterior probability can be 
re-written as 



P{F = x,\N) 



PjNlF ^ X,) 

160 

J2 P{N\F = x^) Ax^ 
1=1 



(3) 



where i runs over the array of simulated flux levels and 
Ax, is step size between each consecutive flux level. 
This posterior probability curve is normalised and can 
be integrated yielding a cumulative posterior probabil- 
ity, P{F < Xui\N)^ from which one can determine the 



flux upper limit, Xui, at a chosen confldence level. Fig- 
ure [9| shows the posterior probability and cumulative 
distribution functions for three different searches. The 
flux value where the cumulative probability distribution 
crosses 0.95, marks the upper limit on the emitted flux 
correlated with GRPs at the 95% confidence level. 

The upper limits on the flux correlated with GRPs at 
the 95% confidence level are plotted in Figure [10] for 64 
of the 72 searches. On the shortest time scales probed, 
the duration of one pulsar period, a limit of ~5-10 times 
the average Crab pulsar flux is set on the inter-pulse and 
the combined inter-pulse and main pulse searches. For 8 
of the searches performed around main pulse GRPs, the 
posterior probability curve had not converged to zero be- 
fore a flux value of 50 times the measured flux. This 
is due to the low rate of main pulse GRPs detected 
at 8.9 GHz, which resulted in a small number of se- 
lected gamma-ray events in some of the main pulse GRP 
searches. The relatively large fluctuations inherent in 
small-number statistics result in wide posterior proba- 
bility curves and, thus, very high upper limit values. We 
do not quote a 95% confldence level upper-limit value for 
these searches, but note that it is likely around 50 times 
the average gamma-ray flux. As the search window size 
is increased, and thus the statistical sample increases, 
these fluctuations decrease and the upper-limit values 
become more constraining. Once enough of the gamma- 
ray sample is selected due to the increasing size of the 
search windows, the upper-limit values level out at ^2- 
2.5 times the average Crab pulsar flux. We note that for 
the longest two search windows, 729 and 2187 periods, 
between 70% and 80% of the VERITAS events are se- 
lected in the interpulse and combined searches resulting 
in a reduced enhancement sensitivity in these searches 

B 

•^^ The 70% and 80% values were determined without inclusion 
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Figure 10. 95% confidence level upper limit on the gamma-ray 
fiux correlated with GRPs for 64 of the 72 searches. We do not 
quote an upper-limit value for 8 searches with short windows placed 
around main pulse GRPs, but note that it is likely around or above 
50 Crab pulsar flux units. The slight negative and positive shift of 
the x-position of the before and after symbols is done as a visual 
aid to prevent clutter about the common x-coordinate value for the 
before, centered and after symbols. 

5.3. Estimation of Uncertainty 

As stated earlier, the fraction of Monte-Carlo searches 
which yield the same number of correlated events as were 
fomrd in the data, P{N\F), was determined from a Pois- 
son fit to the distributions, rather than from the dis- 
tributions themselves. By adopting this procedure we 
found that fluctuations in the posterior probability den- 
sity curves are dramatically reduced, while having little 
effect on the computed 95% confidence level upper-limit 
values. We investigated the uncertainty on the poste- 
rior probability values due to the finite statistics in the 
Monte-Carlo distributions and the Poisson fitting proce- 
dure. 

The values were determined between a given Monte- 
Carlo distribution and Poisson curves with a range of 
mean values centered on the best fit Poisson mean value. 
From these values, a fit probability versus Poisson 
mean curve was generated, and from this curve, random 
values were drawn. A histogram of P{N\F) values was 
then compiled using Poisson curves with these random 
values as their mean. The standard deviation of this his- 
togram was then used as the uncertainty on the P{N\F) 
value determined from the best fit Poisson curve. The 
uncertainty on the probability values determined in this 
way is less than 15% for the bulk of distributions but is 
occasionally as large as 30%. The error bars on the prob- 
ability values plotted in Figure |9] were calculated in this 
way. Other methods to calculate the uncertainty were 
investigated and found, in general, to yield a smaller un- 
certainty value. Folding the uncertainty on the probabil- 
ity values into the cumulative distribution, however, the 
uncertainly on the 95% confidence level upper-limit value 
was found to be less than 3% for the 64 searches which 
yielded a limiting value. This underlines the robustness 
of the procedure we used to calculate the upper-limit 
values. 

of the data acquired on the first day of observations where our 
GRP count was relatively low. 



5.4. Selecting only the most energetic GRP events 

Having investigated GRPs whose peak flux density is 
greater than 7a above the averaged radio signal and ob- 
served no VHE enhancement, we now consider only the 
most energetic GRP events. Due to the low rate of GRP 
events in the 2008 data set, the following analysis uses 
the data collected in 2009 only. 

Earlier, we estimated that the system equivalent flux 
density for the 2009 observations was 0.7 Jy. The sam- 
pling interval for these observations was 3.2 fj,s. Thus, 
the energy of a GRP is 

GRPe = ^^^(0.7Jy) X (iV,an,p)(3.2//s) (4) 

where Gi?Pg/N is the signal-to- noise ratio of the GRP 
and -/Vsamp i s the width of the GRP in samples (see 
IBilous et all (I20TI) . for more details regarding the en- 
ergy calculation of GRPs). The effective width of each 
GRP is determined from the width of the boxcar function 
which, when convolved with the radio time series, returns 
the largest S/N value for a given pulse. The widest GRP 
found had a width of 14 samples, or 44.8 ^s. 

In order to have a common energy scale for observa- 
tions taken on different nights, we must correct for the 
effects of RISS. This was done by compiling average ra- 
dio pulse profiles on 80-minute time scales, which is the 
characteristic time scale of RISS for the Crab at 8.9 GHz. 
The profile for the second half of the observing session 
on MJD 55180 was the strongest and we chose it as ref- 
erence pulse profile to which the other profiles were com- 
pared. Ratios of the peak S/N values were computed 
between the reference profile and the other profiles and 
used as RISS-correction coefficients to scale the GRP en- 
ergy values computed by Equation [?1 Figure [TT] shows 
the resulting distributions of GRP energy, which follow 
the expected power-law shape with a spectral index of 
—4.03. The roll-off below ~60 Jy/zs is caused by a bias 
in our GRP selection, which was based on peak fluxes 
and reduces the number of broad weak pulses selected. 
The different number of GRPs selected below the roll- 
off on each night of observation is due to the different 
amount of RISS contributions to the GRP fiux densities 
on a given night. 

Having computed an RISS-corrected energy value for 
every GRP event in the 2009 data set, we repeated the 
correlation search with the VERITAS data. We applied 
three different energy cuts to the radio events, namely: 
60 Jy^s, 100 Jy/xs, and 150 Jy/zs. The 60 Jy^s thresh- 
old was chosen as a low common energy threshold which 
provides an unbiased sample of GRPs, independent of 
both the GRP search bias and the RISS conditions dur- 
ing a particular observing session. The other two thresh- 
olds were chosen to be "high" and "very high" , selecting 
the most energetic ^6% and '-^1% of GRPs, respectively. 
The results of these searches are plotted in Figure [T2l 
We find no significant enhancement (> 3(t) in the VHE 
emission correlated with these energetic GRP events. 

6. DISCUSSION AND CONCLUSION 

Following our study of simultaneous radio and gamma- 
ray data we observe no significant enhancement in VHE 
gamma-ray emission from the Crab pulsar correlated 
with GRPs observed at 8.9 GHz. Our findings are sim- 
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Figure 11. The distribution of giant pulse energies for the latter 
three nights of observations. The distributions follow a power law, 
with the combined data set having a spectral index of —4.03. The 
different number of GRPs selected below the ~60 Jy/^s roll-off on 
each night of observation is due to the different amount of RfSS 
contributions to the GRP flux densities on a given night. 

ilar to those previo usly reported at lowe r Kamma-ray 
and X-ray energies (jLundgren et al.l 119951 : iBilous et al.l 
I2011L [2013 1. Given the level of uncertainty in theories of 
GRP emission, it is hard to draw firm conclusions result- 
ing from the lack of any observed correlation with VHE 
emission. We are not aware of any theory with quanti- 
tative predictions of correlated emission between GRPs 
and VHE emission. GRP-emission mechanisms associ- 
ated with changes in plasma coherence will not cause en- 
hancements in incoherent emission. Small and localised 
changes in the pair-creation rate, which can explain the 
sm all (3%) optical enh ancements previously measured 
bv lShearer et all (|2003| ). would yield VHE flux enhance- 
ments which are below our sensitivity. 

Enhanced gamma-ray emission in connection with 
high-frequency interpulse GRPs, as postulated by 
iLvutikovl ()2007[) . is not observed. In this model the 
primary beam, which generates coherent GRPs through 
anomalous cyclotron resonance, also emits curvature 
gamma rays with energies in the tens of GeV range. 
It has since been realised that there is a solid upper 
limit on the en ergy of curvature rad iation in the Crab, 
< 150 GeV (|Lvutikov et al.llMl) . This is due to the 
equivalence of the rate of energy loss and rate of accel- 
eration gains achieved by charged particles energised in 
the outer magnctosphcrc. Thus, in the Crab pulsar, the 
curvature photons from the primary beam generally do 
not reach the VHE energy band. This means that the 
non-detection of a VHE enhancement, presented here, 
does not contradict this GRP emission model. 

However, it is possible that there is an indirect link 
between VHE and radio emission for interpulse GRPs. 
Recent models for the pulsed VHE gamma-rays detected 
from the Crab suggest that inverse- Compton emission 
dominates curvature emission at energies above a few 
GeV in the ou ter magnetosphere (jLvutikov et al.l [20121 : 
iDu et al.1l2012l) . A recent study of the Geminga pulsar 



also supports this scenario ()Lvutikovl[20T2l ). In these 
models, the primary particle beam is accelerated in a 
modest electric field (of few percent of the magnetic field 
strength) in the outer magnetosphere and produces cur- 
vature emission up to '^10 GeV. Photons in the VHE 
band are generated by secondary pairs which up-scatter 
their own cyclotron or synchrotron emission. Within the 
framework of the high- f reque ncy interpulse GRP emis- 
sion model of ILvutikovl (|2007[ ). a VHE-GRP connection 
is thus still expected. In this model, roughly half of 
the energy of the plasma beam, energised in the mag- 
netic reconnection event, goes into the production of a 
secondary plasma which can in turn generate enhanced 
VHE gamma-ray emission via synchrotron self- Compton 
scattering. The level of the VHE enhancement is, how- 
ever, difficult to estimate given the uncertainties and 
fluctuations in the properties of both the primary and 
secondary beams. This enhancement is likely below the 
sensitivity of present VHE instruments. 

In the model of pulsed V H E em ission from the Crab 
pulsar of lAharonian et al.l (|2012[ ). where pulsed X- 
ray emission originating in the magnetosphere is up- 
scattered to VHE energies in the wind beyond the light 
cylinder, a GRP- VHE correlation should exist at most 
at the same level as a GRP-X-r ay correlation. Sin ce no 
GRP-X-ray correlation is seen (jBilous et al.l 120121 ) and 
the upper limit on the enhanced X-ray flux is 30%, any 
corresponding VHE enhancement is expected to be below 
30% and thus below our sensitivity. 

This research is supported by grants from the U.S. 
Department of Energy Office of Science, the U.S. Na- 
tional Science Foundation and the Smithsonian Institu- 
tion, by NSERC in Canada, by Science Foundation Ire- 
land (SFI 10/RFP/AST2748) and by STFC in the U.K. 
Wc acknowledge the excellent work of the technical sup- 
port staff at the Fred Lawrence Whipple Observatory 
and at the collaborating institutions in the construction 
and operation of the instrument. The National Radio 
Astronomy Observatory is a facility of the National Sci- 
ence Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

REFERENCES 



Abdo, A. A., et al. 2010, ApJ, 708, 1254 

Aharonian, F. A., Bogovalov, S. V., & Khangulyan, D. 2012, 

Nature, 482, 507 
Aleksic, J., et al. 2012, A&A, 540, A69 
Aliu, E., et al. 2011, Science, 334, 69 

Argyle, E., Baird, G., Grindlay, J., Ifelmken, If., & Omongain, E. 

1974, Nuovo Cimento B Serie, 24, 153 
Bilous, A. v., Kondratiev, V. I., McLaughlin, M. A., Ransom, 

S. M., Lyutikov, M., Mickaliger, M., & Langston, G. I. 2011, 

ApJ, 728, 110 

Bilous, A. v., McLaughlin, M. A., Kondratiev, V. I., & Ransom, 

S. M. 2012, ApJ, 749, 24 
Burke-Spolaor, S., et al. 2012, MNRAS, 423, 1351 
Cordes, J. M., Bhat, N. D. R., Hankins, T. H., McLaughlin, 

M. A., & Kern, J. 2004, ApJ, 612, 375 
Davies, J. M., & Cotton, E. S. 1957, Solar Energy, 1:2-3, 16 
Du, Y. J., Qiao, G. J., & Wang, W. 2012, ApJ, 748, 84 
Fierro, J. M., Michelson, P. F., Nolan, P. L., & Thompson, D. J. 

1998, ApJ, 494, 734 
Fomin, V. P., Stepanian, A. A., Lamb, R. C, Lewis, D. A., 

Punch, M., & Weekes, T. C. 1994, Astroparticle Physics, 2, 137 
Hankins, T. H., & Eilek, J. A. 2007, ApJ, 670, 693 



12 



10 r 



10-^ 



10 - 



: (a) - Main pulse > 60Jy[iS 














i *^ 












■ " : : 'i! " ; 


A Befo 






■ Centered 


• A ■ 
■ • -■- 


• After 
— MC Range 



i 



iir "I 



(b) - Interpulse > 60Jyp.s 



A 


Before 


■ 


Centered 


• 


After 




MC Range 



H+t| 1 — I I I I I ll| 1 — I I I llll| 1 — I I I 1 1 ll| 1 — h 

(e) - Interpulse > lOOJyus 



(c) - Combined > 60Jyns 



A 


Before 


■ 


Centered 


• 


After 




MC Range 



t+tt] 1 — I I I lllll 1 — I I I lllll 1 — I I I lllll 1 — I- 



« 10' - 



10 r 



1 =- 



; (d) 


- Main pulse > lOOJyus 








V 










J .1. II 


1 ^ 1 


A Before 
■ Centered 
• After 
— MC Range 









• 



A Before 
■ Centered 
• After 
— MC Range 



: (f) - Combined >100Jy|iS 












of 










ii i 




A Before 
■ Centered 

• After 
— MC Range 







" I' T T I 1 1 I I II 

(g) - Main pulse > ISOJyus 



t+ttf- — I t I 1 1 lll| 1 I I I llll| 1 I I I llll| 1— h 

(h) - Interpulse > ISOJyus 



Htif- — I I I I llll| 1 I I I llll| 1 I I I llll| 1— I- 

(i) - Combined > ISOJyus 



10^ r 



lO'' 



10 



III III fli 



, Before 
I Centered 
< After 
- MC Range 



III III 







• 




A Before 
■ Centered 
• After 
— MC Range 



III 



A 


Before 


■ 


Centered 


• 


After 




MC Range 



I 



I 



I 



I I I 



1 10 10 10 

Search Window Duration [Pulsar Periods] 



1 10 lO-" 10 

Search Window Duration [Pulsar Periods] 



1 10 10^ 10 

Search Window Duration [Pulsar Periods] 



Figure 12. The number of gamma-ray events selected (symbols) versus the duration of the search window. The grey regions denote 68% 
containment intervals about the mean of the distribution (black line) determined from the identical searches performed on Monte-Carlo 
data sets containing an injected signal from the Crab pulsar at the level of the measured VHE pulsar flux. The upper, middle and lower 
rows show the results of the enhancement search when restricting the energy of the GRP events to be above 60 Jy/^s, 100 Jy/^s and 150 Jy^ts, 
respectively. No excess (or deficit) is found in any search with a probability equivalent of 3a or higher. See the caption of Figure [6] for a 
further description of this figure. 
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